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Abstract 

We discuss charge symmetry and charge independence breaking in a chiral effective field theory 
approach for few-nucleon systems based on a modified Weinberg power counting. We construct a 
two-nucleon potential with bound and scattering states generated by means of a properly regu- 
larized Lippmann-Schwinger equation. We systematically introduce strong isospin-violating and 
electromagnetic operators in the theory. We use standard procedures to treat the Coulomb poten- 
tial between two protons in momentum space. We present results for phase shifts in the proton- 
proton, neutron-proton and the neutron-neutron systems. We discuss the various contributions 
to charge dependence and charge symmetry breaking observed in the nucleon-nucleon scattering 
lengths. 
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1 Introduction 



The forces between two (or more) nucleons have been studied over many decades but the connection to 
the underlying theory of the strong interactions is so far only loosely established. QCD is formulated 
in terms of strongly interacting quark and gluon fields and thus not amenable to most standard field 
theoretical methods.^ On the other hand, it is well established that meson theory works fairly 



well at low energies. Such models are based on the exchange of pions and heavier mesons between 
the nucleons supplemented by form factors or short distance regulators. The drawback with such 
approaches is the fact that one can not systematically calculate corrections and that strong form 
factors are ill-defined concepts in any field-theoretical sense. However, the pion exchange originally 
proposed by Yukawa can now be linked to one specific property of low-energy QCD, namely its 
spontaneous and explicit chiral symmetry violation. Pions arise as Goldstone bosons with point- 
like derivative couplings to the nucleons in accordance with Goldstones theorem. In fact, one has a 
hierarchy of scales in the problem, the long distances given by one- and two-pion exchange, whereas 
intermediate and short distances are often modelled in terms of heavy meson exchanges and also 
are given by the unphysical form factor ranges. Starting from nucleoli and pion degrees of freedom, 
effective field theory (EFT) allows for a clean separation of these scales and consequently offers a 
systematic and controlled understanding of the forces between two or more nucleons. Apart from 
dealing with the various scales appearing in nuclear systems, it is straightforward to implement the 
spontaneously and explicitely broken chiral symmetry of QCD as well as external probes in the EFT. 
The appearance of shallow bound states (or, equivalently, large scattering lengths) requires some 
method of resummation. There are essentially two ways of tackling this problem. One approach is 
to build the potential from EFT and employ it in a properly regularized Lippmann-Schwinger (or 
Schrodinger) equation. This scheme was proposed in ref.[|J and will be referred to as the Weinberg 
approach. Alternatively, one can also do the expansion directly on the level of the scattering amplitude 
and resum the leading order, momentum-independent four-nucleon interaction This is the so- 
called KSW scheme. While at extremely low energies, much below the scale set by the pion mass, it 
is sufficient to consider four-nucleon interactions only (otherwise the rather successfull effective range 
expansion would not work), for typical nuclear momenta of the size of the pion mass, pions have to 
be included explicitely. In refs.[||, |J a next-to-next-to-leading order chiral two-nucleon potential 
in a modified Weinberg scheme was developed for the neutron-proton {np) system and shown to be 
close in accuracy to the so-called modern potentials (in some partial waves). This isospin symmetric 
potential will be the starting point of the investigation presented here. 

It is well established that the nucleon-nucleon interactions are charge dependent (for a review, see 
e.g.|| and a recent summary is given in ||). For example, in the 1 So channel one has for the scattering 
lengths a and the effective ranges r 

a C iB = - (a nn + a pp ) - a np = 5.64 ± 0.40 frn , (1.1) 
rem = ^ (r n „ + r pp ) - r np = 0.03 ± 0.06 fm . (1.2) 

These numbers for charge independence breaking (CIB) are taken from the recent compilation of 
Machleidt (?]]. It is understood that the Coulomb effects for pp (nn) scattering are subtracted based 
on standard methods (for a treatment of the Coulomb effects in proton-proton scattering in an 
EFT framework, see ref.[H). One sometimes refers to these as the nuclear amplitudes. The charge 

* 4 In principle, lattice gauge theory can be used to study this problem but so far it has not attracted much attention. 
In any case, a purely numerical solution to the problem always needs to be supplemented by additional physical insight 
obtained using other methods. 
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independence breaking in the scattering lengths is large, of the order of 25%, since a np = (—23.740 ± 
0.020) fm. Of course, it is magnified at threshold due to kinematic factors (as witnessed by the 
disappearance of the effect in the effective range)^. In addition, there are charge symmetry breaking 
(CSB) effects leading to different values for the pp and ran threshold parameters, 

ocsb = a pp - a nn = 1.6 ± 0.6 fm , rcsB = r pp - r nn = 0.10 ± 0.12 fm . (1.3) 

Combining these numbers gives as central values a nn = —18.9 fm and a pp = —17.3 fm. We should point 
out already here that a recent measurement performed at Bonn |J leads to a sizeably smaller value for 
a nn in the 1 Sq channel which affects both values ocib and acsB- We come back to the implications of 
that measurement later. Both the CIB and CSB effects have been studied intensively within potential 
models of the nucleon-nucleon (NN) interactions. In such approaches, the dominant CIB comes from 
the charged to neutral pion mass difference in the one-pion exchange (OPE), a^f^ ~ 3.6 ± 0.2 fm. 
Additional contributions come from 777 and 2ir (TPE) exchanges. According to some calculations, 
their size is approximately l/3 rd of the leading OPE contribution. For the TPE, significantly smaller 
results can also be found in the literature, see e.g. table 3.3 in ref.(||. In case of the ttj contribution, 
a recent calculation[10] gives a value which is more than a factor of ten smaller than the leading 



OPE contribution. We will come back to these issues later. Note also that the charge dependence in 
the pion-nucleon coupling constants (if existing) in OPE and TPE almost entirely cancel. Naively, 
one would expect a possible charge dependence of the pion-nucleon coupling constants to be of a 
similar importance as the pion mass difference in OPE. This strong suppression of charge-dependent 
couplings has eluded a deeper understanding until the EFT considerations based on power counting 



spelled out in ref.[ll]. In the meson-exchange picture, CSB originates mostly from p — uj mixing, 
a CSB ~ 1-2 i 0.4 fm. Other contributions due to tt — 77, vr — 77' mixing or the proton-neutron mass 



difference are known to be much smaller. This picture has been challenged in ref.[12] where it was 
claimed that the sum of various two-boson-exchange diagrams can fully explain the observed CSB 
in the singlet scattering length (for a more detailed discussion and a comparison to the approaches 
based on p — to mixing, see ref . ) . It was already shown in refe.0, 0, ^ that EFT 

gives novel 

insight concerning the size of the TPE and 7T7 contributions as well as the suppression of possible 
charge-dependent irN coupling constants. However, a systematic and detailed study of CIB and CSB 
in the Weinberg approach has not been given so far. We will fill this gap by working out the complete 
effects of isospin violation at leading order (LO) and next-to-leading order (NLO) in the modified 
Weinberg approach of refs.||, 0], which has been proven very successful up to the pion production 
threshold, i.e. E Uh ~ 300 MeV. 

The manuscript is organized as follows. In section || we give a brief reminder of the effective field 
theory underlying our calculations. The problem of isospin violation in the NN force is taken up 
in section ||| We show how to incorporate strong and electromagnetic isospin breaking in the EFT 
based on an extended power counting and organize the various contributions from the light quark 
mass difference as well as from soft and hard photons. In addition, we discuss how to treat the 
Coulomb potential in momentum space. In section |] the renormalized potential at NLO is displayed. 
It consists of one- and two-pion exchanges, 717 graphs and local four-nucleon interactions with zero 
or two derivatives. The results are given and discussed in section ||. We conclude with a summary 
and outlook. Some subtleties related to the CIB two-pion exchange are spelled out in the appendix. 



* J Sometimes, CIB is simply discussed as the difference between the Coulomb subtracted pp and the np observables, 
which has to be accounted for when comparing to other values obtained in the literature. To avoid confusion, we will 
always explicitely state the definition we are using 
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2 Brief reminder of nuclear chiral perturbation theory 



In this section we briefly spell out the central ideas underlying our calculations. One starts from an ef- 
fective chiral Lagrangian of pions and nucleons, including in particular local four-nucleon interactions 
which describe the short range part of the nuclear force, symbolically 

£eff = -Ctttt + C-ttN + C-NN , (2.1) 

where each of the terms admits an expansion in small momenta and quark (meson) masses. To a 
given order, one has to include all terms consistent with chiral symmetry, parity, charge conjugation 
and so on. The methods to construct such Lagrangians and to calculate the corresponding Feynman 



diagrams are by now standard, see e.g. [16]. From the effective Lagrangian, one derives the two- 
nucleon potential. This is based on a modified Weinberg counting as described in ref.Q which is 
applied to the two-nucleon potential to a certain order in small momenta and pion masses, 

V(p,p')=J2V (i) (p,p') , (2-2) 

i 

with p, p' the nucleon centre-of-mass (cms) momenta and the superscript i gives the (non-negative) 
chiral dimension. The power counting underlying this potential is based on the considerations pre- 
sented in ref.||]. To leading order (LO), this potential is the sum of the time-honored one-pion 
exchange (OPE) (with point-like coupling) and of two four-nucleon contact interactions without 
derivatives. The low-energy constants (LECs) accompanying these terms have to be determined by a 
fit to some data, like e.g. the two S-waves in the low energy region (for np). At next -to-leading order 
(NLO), one has corrections to the OPE, the leading order two-pion exchange (TPE) graphs (called 
the box, triangle and football diagrams according to the number of non-linear NNtttt interactions 
involved, which are zero, one and two, in order) and seven dimension two four-nucleon terms with 
unknown LECs (for the np system). The existence of shallow nuclear bound states (or large scatter- 
ing lengths a 1/M n ) forces one to perform an additional nonperturbative resummation. This is 
done here by obtaining the bound and scattering states by solving a regularized Lippmann-Schwinger 
(LS) equation, projected onto states with orbital angular momentum I, total spin s and total angular 
momentum j, 

Jo yiTi) ' p A — p" +%7] 

with n — > + and m denotes the nucleon mass. In the uncoupled case, I is conserved. The partial 
wave projected potential V^ip' , p) is obtained by standard methods, see e.g. [|J|. The potential has 
to be understood as regularized, as dictated by the EFT approach employed here. That is done in 
the following way: 

V(p,p' ) - f R (p) V(p,p') f R (p') , (2.4) 

where fii(p) is a regulator function chosen in harmony with the underlying symmetries. In what 
follows, we work with a sharp cut-off, 

f^(p) = 9(A 2 -p 2 ). (2.5) 

Within a certain range of cut-off values, the physics should be independent of its precise value. This 
range increases has one goes to higher orders, as demonstrated explicitely for the np case in H. 
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3 Explicit isospin violation 



3.1 Sources of isospin violation 

Here, we briefly discuss how isospin violation arises in the Standard Model and show how it can be 
implemented in the effective field theory. We stress that there are two small parameters which allows 
one to treat strong and electromagnetic isospin breaking perturbatively. The electromagnetic effects 
can further be separated into two types, soft and hard photons. The soft photons lead on one hand to 
perturbatively calculable photon loop effects and on the other hand the nonperturbative resummation 
of such photons generates the long range Coulomb force. The effects of hard photons can completely 
be absorbed in electromagnetic short distance operators. We also briefly discuss the relation to the 
notions of charge independence and charge symmetry more often used in this context. 

Consider first pure QCD, i.e. a world with no electroweak interactions.^ To be specific, we are 



only concerned with the sector of the two light up and down quarks. For equal current quark masses 
m u = md, the QCD Hamiltonian is invariant under a global (flavor) transformation of the type 



Vq = V " , VeSU(2) . (3.1) 





Of course, in nature the light quark masses are not equal. In fact, they differ considerably, m u /m,i = 
0.55 ±0.15 for the standard MS subtraction scheme at a renormalization scale of 1 GeV. Still, isospin 
is a good approximate symmetry because what counts is not the ratio of the quark masses but rather 
their difference compared to the typical scale of strong interactions, 

£ _ d «~ «!, (3.2) 

Aqcd 30 

since Aqcd ~ 150 MeV. QCD thus provides a reason why isospin is such a good approximate sym- 
metry. In fact, one could also use the larger chiral symmetry breaking scale A x = 4TrF n / \fNf, with 
F n = 92.4 MeV the pion decay constant and Nf the number of quark flavors. This would indicate even 
smaller isospin breaking effects. For our purposes it suffices to state that such dimensional arguments 
let one understand the smallness of the observed isospin breaking. Another way of looking at the 
relative strength of isospin violation is to compare the isoscalar (isospin conserving) ~ (m u + m^) to 
the isovector (isospin violating) ~ (m u — rrid) part of the QCD quark mass term, 

W2S = ~?(mu + m d )(l + 673)9, (3.3) 

leading one to expect that 

m d -m u 1 

e = ■ 3.4) 

m d + m u 3 

would be the appropriate small parameter. From the numerical value of e, isospin symmetry would 
appear to be accidental. This argument, however, ignores the generic scale of QCD generated by 
dimensional transmutation. On the other hand, for the organisation of the effective field theory, it 
is advantageous to work in the isospin basis as given in eq.( |3.4| ). Therefore, we will count strong 
isospin violation in terms of e, but keep in mind that the true size of isospin violation is smaller than 
indicated by the numerical size of the parameter e. 

In the presence of electromagnetism, further isospin violation is induced since the charges of the 
quarks are unequal. These effects are in general small due to the explicit appearance of the fine 



#6 We consider the quark masses as parameters and are not concerned with their dynamical origin here. 
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structure constant a = e 2 /47r = 1/137.036. In fact, in many cases the strong and electromagnetic 
isospin-violating effects are of the same size. A good example is the neutron-proton mass difference 
where (m n — m p )^ CD ~ 2.1 MeV and (m n — m p )^ ED ~ — 0.8MeV, leading to the physical value of 
m n — m p = 1.3 MeV. The pion mass difference is special, because of the absence of D-couplings in 
SU(2) there is no strong contribution in the light quark mass difference but only a very tiny second 
order contribution (due to 7r° — 7] mixing). In what follows, we will use 

AM^ = M n ± - M w0 = 4.6 MeV , (3.5) 

which is almost entirely of electromagnetic origin. Due to the long range nature of the electromagnetic 
interactions, one has to be careful in separating long and short distance physics generated by virtual 
photons.^ 7 While the latter can be represented by a string of local contact interactions with increasing 



chiral dimension, the former generate the Coulomb singularities and need to be treated separately. 

In the nuclear physics language, it is more common to talk about charge (in) dependence and charge 
symmetry (breaking). Charge independence refers to the invariance under any rotation in isospin 
space. A violation of this symmetry is called charge dependence or charge independence breaking 
(CIB). On the quark level, charge symmetry is a rotation about the 2-axis in isospin space 

"-(;)=^(2)=(t) • • ^ 

and similarly for the nucleon isodoublet consisting of the proton and the neutron. The violation of 
this symmetry is called charge symmetry breaking (CSB). It is a special case of charge dependence. 
On the level of the two-nucleon force, charge independence implies a potential of the form (i.e. which 
is a scalar in isospin space) 

Vij =A + Bn- tj , (3.7) 

where A and B are functions of the nucleon spin and momentum (or space) coordinates and 
labels the interacting nucleons. Charge symmetry allows for a more general two-body force, 

V ij)CS =A + Bn- Tj + Crf ■ rf . (3.8) 

If one works out the Coulomb potential in terms of the nucleon charge matrix Q, it is obvious that 
the electromagnetic effects lead to breaking of charge independence and charge symmetry, 

e 2 

^Coulomb = t^— 2 (1 + rf)(l + r 3 ) . (3.9) 

Stated differently, CIB of the strong NN interactions means that the proton-proton (pp), neutron- 
proton (np) or neutron-neutron (ran) are different once the Coulomb effects have been removed. CSB, 
on the other hand, refers to the difference between the pp and the nn interactions. 



3.2 Power counting, effective Lagrangian and classification scheme 

As discussed in the preceeding section, the Standard Model supplies us with two small parameters, 
which allow one to systematically include isospin violation in terms of two external scalar sources, 

X = 2B (m u + m d ) (1 + er 3 ) , (3.10) 
Q = f(l + r 3 ), (3.11) 
* 7 By virtual photons we refer to photon exchanges within any given Feynman diagram. 
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where B measures the strength of the quark-antiquark condensation in the vacuum, B = |(0|uu|0)|/-F 2 , 
with F n the weak pion decay constant. We work here in the standard approach with B ^> F n . The 
charge matrix Q defines the strength of the virtual photon coupled to the matter fields. Since a 
virtual photon can not leave a Feynman diagram, only even powers of Q can appear in the effective 
Lagrangian (with the exception of the covariant derivative), as detailed below. A further consequence 
is that virtual photons appear necessarily in loop diagrams. The germane electromagnetic expansion 
parameters is a = e 2 /4ir (although one could equally well use e 2 or e 2 /(4-7r) 2 ). 

Due to its perturbative nature induced by the small parameters, we treat the strong and electromag- 
netic isospin violation in addition to the power counting of the isospin symmetric potential mentioned 
in section ||. In principle, one has a double expansion of the form e n a m with n, m are integers. How- 
ever, since we iterate the potential in the LS equation, we have to avoid double counting. Here, the 
isospin basis given in Eqs.( |3.10 , |3.1lD proves to be very useful. 



We now discuss the various parts of the effective Lagrangian underlying the analysis of isospin violation 
in the two-nucleon system. To include virtual photons in the pion and the pion-nucleon system 
using the above-defined external sources is by now a standard procedure |17]-[24|. The lowest order 
(dimension two) pion Lagrangian takes the form 

Anr = tfiVpUVrt + X U ] + x ] U) + C(QUQU^) . (3.12) 

Here, U(x) collects the isotriplet of pion fields, is the (pion) covariant derivative containing 
the virtual photons, ( ) denotes the trace in flavor space, and the last term, which contains the 
nucleon charge matrix Q=ediag(l, 0) J^l leads to the charged to neutral pion mass difference, AM 2 = 
M 2 ± - M 2 , via AM 2 = 8iraC/F%, i.e. C = 5.9 • 1(T 5 GeV 4 . Note that to this order the quark mass 
difference m u — does not appear in the meson Lagrangian (due to G-parity). That is chiefly the 
reason why the pion mass difference is almost entirely an electromagnetic (em) effect. The equivalent 
pion-nucleon Lagrangian to second order takes the form 




1 2m I 8m 



+ \ ( C4 + zL") ^ a ^ UiU i + C5 (*+ " £<x+>) + •••} N ' ( 3 - 13 ) 

which is the standard heavy baryon effective Lagrangian in the rest-frame = (1,0,0,0). m is the 
nucleon mass, u = y/U, the chiral viel-bein, ~ —id^/F^ + . . . and \+ = ux'u + v)x u ^ ■ The 
four-nucleon interactions to be discussed below do not modify the form of this Lagrangian (for a gen- 
eral discussion, see e.g. ref.[22|). Strong isospin breaking is due to the operator ~ eg. More precisely, 



this low-energy constant always appears in the combination C5 B (m u — m^). Electromagnetic terms 
to second order are given by [23, 24] 



-nN 



F 2 iVt [f^Ql - Q 2 _) + f 2 Q + (Q + ) + MQl + Ql)} N , (3.14) 
with Q± = uQ^u^u^Qu^ and A = A — (A)/2 projects onto the off-diagonal elements of the operator 



A. The last term in eq.(3.14) is not observable since it leads to an equal electromagnetic mass shift 
for the proton and the neutron, whereas the operator ~ f 2 to this order gives the em proton-neutron 
mass difference. In what follows, we will refrain from writing down such types of operators which only 
lead to an overall shift of masses or coupling constants. We note that in the pion and pion-nucleon 



* 8 Or equivalently, one can use the quark charge matrix e(l + T3)/2, see refs.|23|, |24|. 
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sector, one can effectively count the electric charge as a small momentum or meson mass. This is 
based on the observation that M 7T /A X ~ e/v47r = \fa ~ 1/10. It is thus possible to turn the dual 
expansion in small momenta/meson masses on one side and in the electric coupling e on the other 
side into an expansion with one generic small parameter. As noted before, we use the fine structure 
constant a = e 2 /Att as the electromagnetic expansion parameter. We now turn to the two-nucleon 
sector, i.e. the four-fermion contact interactions without pion fields. Consider first the strong terms. 
Up to one derivative, the effective Lagrangian takes the form 

£f N = h(NiN) 2 + l 2 (N^aN) 2 + l 3 (N^( X+ )N)(N^N) + h(rt x+N)(rf N) 

+ l 5 (N^a( X+ )N){N^aN) +k(N^ax+N)(N^aN) + . . . , (3.15) 

where the ellipsis denotes terms with two (or more) derivatives acting on the nucleon fields. Similarly, 
one can construct the em terms. The ones without derivatives on the nucleon fields read 

{ri(Ql ~ Q-) + r 2 Q+{Q+}} N(N^N) 

N^a {r 3 (Q 2 + - Q 2 _) + r 4 Q + (Q+)} N(N^aN) 

+ Aft { r5 Q + + re (Q + )} N(N^Q + N) + rfa{r 7 Q+ + r 8 (Q+)} N(N^aQ + N) 

+ r 9 (N^Q + N) 2 + r 10 (N^aQ + N) 2 . (3.16) 

There are also various terms resulting form the insertion of the Pauli isospin matrices f in different 
N^N binomials. Some of these can be eliminated by Fierz reordering, while the others are of no 
importance for our considerations. The special form of eqs. (3. 10 j3.ll ) allows us to rewrite all nucleonic 



-NN — 
+ 



contact terms with insertions of external fields at a given chiral order as the sum of isospin symmetric 
contact structures and of 

Vcsb ~ (N^t 3 N)(N^N), (3.17) 
V CIB ~ (NW 3 N) 2 , (3.18) 

which parameterize the non-pionic CSB and CIB-effects (as discussed in more detail below). Note 
that in the iteration all isoscalar pieces in the external sources Q and x get resummed and we 
therefore only consider terms linear in a and e. As in ref.pl], we form appropriate combinations of 
these four-nucleon operators so that we can pin down the corresponding LECs in the low partial 
waves. Specifically, we combine the isospin conserving and violating terms in the following manner 
(in a somewhat symbolic notation), 

CNN = Cf m (ATt p p N f + C str (Ar t p {x+) N) {N i N) + C orn ^ ^ (q2) n) ^ n) + (3^9) 

where the ellipsis stands for terms with two (or more) derivatives. Also, Pp is a projector onto the 
appropriate partial wave with (3 denoting the corresponding angular momentum and isospin quantum 
numbers, i.e. for (3 = 1 Sq we have Pp = git^t '/V~8- The precise number of parameters in each 
partial wave to a given order will be enumerated in the next section. It is also worth to stress that 
the isospin-breaking contact interactions appear at the same order as the four-nucleon contact terms 
with two derivatives in the isospin symmetric case. For that, we have to specify more explicitely how 
to organize the various isospin breaking contributions to the potential. 

After having set up the power counting of the underlying effective Lagrangian, we are in the position 
of classifying the various contributions to isospin violation (or CIB and CSB) in the NN interaction, 
more precisely, on the level of the two-nucleon potential. The pion mass difference in the OPEP leads 
to a contribution of the form 

V^r ~ AM^Vope ~ ^Vope - aQ~ 2 , (3.20) 
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where Q denotes some small momentum or pion mass and Vope the isospin symmetric OPEP. There- 
fore, the pion mass squared in the denominator leads to an enhancement of this contribution by two 
powers in small momenta. So the leading order isospin violating potential scales as a/Q 2 . This order 
we denote by "L0', meaning leading order in the isospin breaking NN potential. It is also worth 
to point out that we count the Coulomb potential as L0 despite the fact that it stems from an 
infinite sum of photon exchanges. Formally, however, the Coulomb potential a/\q — q'\ 2 has the same 
structure as the L0 OPEP. The leading order momentum-independent isospin violating four-nucleon 
contact interactions do not appear at L0 because they simply scale as aQ° and eQ°, so they come 
together with the pion mass difference and irj graphs at NL0 , since these latter two also scale as 
aQ°. Subleading TPE proportional to the dimension two LECs Cj would appear at N 2 L0 whereas 
a possible isospin breaking in the pion-nucleon coupling constants, g n Opp / g^+pn, is a N 3 L0 effect, 
i.e. is expected to be very suppressed. 



Order 


Parameter 


Contribution 


L0 


a 


Pion mass difference in OPE 








a 


Coulomb potential 






NL0 


a 


Pion mass difference in TPE 








a 


7T7— exchange 








a 


four-nucleon contact interaction with no 


derivatives 


~ (nW 3 n) 2 




€ 


four-nucleon contact interaction with no 


derivatives 


~ (N^t 3 N)(N^N) 



We remark that this power counting is similar to but also distinctively different from the one in 



the KSW scheme [11]. In that approach, pion exchange is treated perturbatively and consequently 
7T7 loop graphs only appear at NNL0. It is instructive to consider the NLO four-nucleon contact 
interactions (which only affect the S-waves). Terms proportional to (N^QN) 2 or (N^x+N) 2 clearly 
lead to CIB since 

(nn\(N^r 3 N) 2 \nn) = {pp\(N^T 3 N) 2 \pp) / (np\ (N^T 3 N) 2 \np) , (3.21) 

whereas the operator proportional to the light quark mass difference, i.e. the isovector part of x leads 
to CSB, 

{nn\(N^T 3 N)(N^ N)\nn) = -(pp|(Ar t T 3 iV)(./V t ./V)|pp) and (np|(iY t r 3 iV)(iV t iV)|np) = . (3.22) 

This means that the leading order mechanisms giving rise to CIB and CSB are very different, in par- 
ticular, CSB is driven by the light quark mass difference. We should also note that in the above classi- 
fication scheme we have not listed the nucleon mass difference explicitely. Throughout, the kinematics 
is given in terms of the physical proton and neutron mass, m p = 938.27 MeV and m n = 939.57 MeV, 
respectively, leading to slightly different threshold energies. This is an effect concerning the external 
legs and therefore affects the S-matrix of the scattering process only indirectly. Consequently, we 
did not list such type of contribution in the classification scheme. In the pertinent TPE graphs we 
neglect this mass splitting, which is consistent with the underlying power counting. We also point 
out that the possible isospin breaking in the pion-nucleon coupling constants appears at NNL0 in 
the KSW scheme as first observed in ref.[|ll]]. Before we can work out the effects of CIB and CSB in 
the two-nucleon system, we have to discuss how to treat the long-range Coulomb force. 



3.3 Momentum space treatment of the Coulomb potential 

As the Coulomb potential is of infinite range, the S-matrix has to be formulated in terms of asymptotic 
Coulomb states. Therefore, the phase shifts for a given angular momentum I due to the strong 
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potential in the presence of the long-range electromagnetic interactions, denoted by <5|, are defined in 
terms of a linear combination of (ir)regular Coulomb-functions F(G) as 

Xi(r) = F l (r)+tan(S l i)Gi(r) (3.23) 

analogously to the expression for an arbitrary potential of short range (i.e. in the absence of the 
Coulomb force) 

Xi(r) = F?(r)+tM8l)G?(r) (3.24) 

with F°, G° denoting solutions of the Coulomb problem with zero charge (conventionally expressed 
in terms of Bessel and Neumann functions) and the corresponding phase shift is called 5f. So far, we 
have restricted ourselves to uncoupled channels, but nevertheless it is straightforward to extend the 
formalism to the case of coupled ones. One simply has to rewrite the wave-functions as 2x2-matrices 

• s = ( V ) • (3 - 25) 

with = F,G, substitute tan<5z by the coupled-channel K-matrix Ki—\ t i+\, and redefine eq.( |3,23| ) as 
a matrix-equation (for compactness, we drop all indices) 

X = F — fiq K G , (3.26) 

with fi being the reduced mass and q the momentum transfer. In what follows we will only consider 
uncoupled channels to simplify notation, nevertheless the extension to coupled channels will always 
be obvious in the above-described sense. 



As eq.( 3,23j ) exhibits asymptotical Coulomb-states, we have to re-express our Lippmann-Schwinger- 



equation in terms of them. Therefore, the calculation of phase shifts requires a large number of 
Fourier-transformed Coulomb-functions. That would make the investigations very time-consuming. 
Consequently, we are now going to develop a formalism that makes use of the phase shifts calculated 
for plane wave asymptotics as defined in eq.( |3.24| ). This procedure has been applied first in ref . 1 25 1 . 
The starting point of this technique is the observation, that for a potential of the form 

V = Vcoulomb + Strong (3.27) 

with 

VkrongV'i = (r > R) , (3.28) 

and vpi the two-nucleon wave-function for a given angular momentum, two exact solutions for the 
wave-function can be given for every point on a sphere with radius R + e. One is of the form as in 
eq.( |3.23| ), and another one according to eq.( |3,24 ) with the phase shifts calculated for the following 



potential as in eq.( 3.27| ), with Vcoulomb being the Fourier-transformed Coulomb-potential integrated 
to the radius R, 



V Cmlomh (\q'-q\) = j R d\e^'-^- 

Jo r 



lKa -(l-co S (\q'-q\R)) . (3.29) 



Here, q, q' are the cms momenta and a is the fine-structure constant. On the above-defined sphere, 
both wave functions describe the same system. Now we know how to obtain an expression for the 
strong shift 5\ in the presence of the Coulomb interaction in terms of the short-range shift 6f in the 
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absence of electromagnetism: We only have to match the two solutions. This is most conveniently 
done by requiring the logarithmic derivative of both solutions to be equal, what enables us to express 
the strong shift in the presence of the Coulomb force in a Wronskian form: 

tan( ^ } " [F ,G] + tan(^)[G ,G] (3 - 30) 

with 

M = (°f-*fL, (3 - 31) 

Here again, the continuation to coupled channels in the form of a matrix-equation is obvious. 

The only remaining difficulty is the determination of the matching radius R, because the given 
solution is wrong as long as ( |3.28| ) is not valid. On the other hand, it is not possible to extend R 
to arbitrarily large values, because the cosine in eq. ( EOsD win cause rapid oscillations. Therefore, 
some test calculations are necessary. First, we put a pure Coulomb potential in our framework. 
In this particular case, the phase shift defined in eq. 

( gjgp should be zero and R should allowed 
to be arbitrarily small (but equal to zero). Nevertheless beyond some radius, the abovementioned 
oscillations are expected. In Fig Jl] we have plotted the phase shift for several lab energies in the 1 5o 
wave. Second, we took the well known Reid93 [^(| potential for proton-proton scattering as V^rong 
and re-did the procedure. In Fig.^ we confer our calculation for different R with the published phase 
shifts of the Reid93-potential at a fixed energy, = lOMeV. As a third check, we compared an 
r-space solution for a Fourier transformed one-pion-exchange potential of the Bonn-type | |27[ | (dipole 
cutoff) and the standard Coulomb potential with the same one-pion-exchange in our framework. The 
checks showed the expected behavior and it seems to be a good choice to define R = 10 fm in all 
partial waves (which is in accordance with a Fourier transformation of our strong potential as well) 
as it has been done for the Bonn potential |]28[| . For R less than 5-7 fm, the strong potential is still 
present, for R bigger than 15 fm, the oscillations can not be neglected. 

3.4 Range expansion 

It is of particular interest to consider the scattering lengths and effective range parameters. The 
effective range expansion takes the form (written here for a genuine partial wave) 

p cot (5) = -- + - rp 2 + v 2 p 4 + v 3 p 6 + v 4P 8 + 0(p 10 ) , (3.32) 
a 2 

with p the nucleon cms momentum, a the scattering length and r the effective range. It has been 



stressed in ref.[29] that the shape parameters Vi are a good testing ground for the range of applicability 
of the underlying EFT since a fit to say the scattering length and the effective range at NLO leads 
to predictions for the higher order range parameters Vi. 

We now turn to the proton-proton scattering length. It is well-known that the strong hadronic 
part a pp can only be separated from the one in the presence of the Coulomb interaction a^j, in a 
regularization scheme dependent way. For a cut-off field theory employing a sharp cut-off as done 
here, one expects a relation of the form, 



1 

+ a 



F ( In — ^) + const. 
m J 



(3.33) 



Some explicit examples for the function in the square brackets in a theory where the pions are 
integrated out are given in ref . ||] (see also [Q ) . We refrain from further discussing the explicit cut- 
off dependence of this relation but will come back to this when we discuss the numerical results. It 
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is also important to stress that within the EFT, only a certain range of cut-off values is allowed and 
thus the often made statement, that the separation of the hadronic scattering length a pp from a pp is 
completely arbitrary, does not hold. 

4 Renormalized potential 

In this section we collect the explicit formulae for the various pieces of the potential expressed in 
terms of renormalized quantities. The infinities appearing in the TPE diagrams are treated exactly 
in the same manner as in ref.[|j] and we refrain from repeating these arguments here. 

Consider first the one-pion exchange (OPE). Its contribution is given in the standard Yukawa form. 
More precisely, for the nn and the pp systems, only the neutral pion can be exchanged whereas for 
the np T = 1 case, we have a superposition of neutral and charged pion exchanges, 

t/(°) _ t/(°) ( 9A \ 2 - - ° l -q^-q n 

^OPE.pp - V OPE,nn ~ ~ [yJT J n ■ T 2 _, 2 + M 2 q ' V 4 --U 

V$kn P = ~ (|r) 2 * ■ r 2 a x ■ qo 2 • <f (^_^ - , (4.2) 

with q = (p' — p) the exchanged momentum in the centre-of-mass system, where p and p ' denote the 
cms momentum of the incoming and the outgoing two-nucleon pair, in order. For later comparison, 
we also give the isospin symmetric one-pion exchange potential, it reads 

t/(°) ( SA \ 2 - -^ox-qo-i-q /a q\ 

Wy m = -( v ^rJ n-Ta , 2 + M , , (4.3) 

with 

M n = -M w+ + -M w o = 138.03 MeV , (4.4) 
3 3 

the average pion mass. So far, we have expressed the pion-nucleon coupling in terms of the Lagrangian 
parameters gA and F n . With the help of the Goldberger-Treiman relation g^NN = 9A m N / F-k it is 
then possible to write down these coupling in a more familiar form. As explained above, we do not 
need to worry about possible isospin breaking effects in g n NN- They are not only suppressed by the 
power counting but also known to be small in the most recent phase shift analysis fll4| . 

Consider now the TPE potential. At NLO, for the case of equal pion masses, we can use eq.(2.12) of 
ref.Q with the average pion mass as defined in eq.( [4.4[) . We give this expression here for completeness, 

= ~3§j^jk L ^ { 4M '( 5 ^ - 4ffi ~ 1) + q\23g\ - 10g\ - 1) + ^ff q2 } 

m = iys?T7i„^±Z±i, (4. 6) 

and we have set q = \q\. Obviously, the first term of the NLO isospin symmetric TPEP in eq.([4.5|) 
is isovector while the second one is isoscalar. The polynom P(k, q) can be absorbed in the isospin 
symmetric dimension two four-nucleon contact interactions up to some small finite shifts as discussed 
in the appendix (in what follows, we will ignore this effect). The pion mass difference in the TPE can 
be incorporated along the lines outlined in ref.|ll|. For completeness, we repeat here the basic steps 
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of that paper and bring the potential in the form appropriate to our discussion. It is most convenient 
to consider the isoscalar and isovector TPE piece separately, 



V2n = vl + V 2 in-T 2 . 

For the isoscalar part, we can express the TPE as 



V, 



2tt 



^Vl{M„ + ,M„ + ) + ±Vl{M^,M^) , 



(4.7) 



(4.8) 



where the arguments refer to the masses of the two exchanged pions. We remark that due to charge 
conservation, only similar pions can be exchanged in this case. This explicit form is used in the 
numerical evaluation with V 2n as given in eq. ( [4.5[) . It is also instructive to perform a Taylor expansion 
of this expression around the average pion mass, 



V&(M„,M„) + 2 



dvl 



2tt 



1 



- <(M 7r ,M 7r ) + 2 



2n 



dM^o 



M 



= M 7r 



(M w + — M n ) + 0(a 2 



--M n 



(M^o - M„) + 0(a 2 ) 



4 dVl 



V&(M )r ,M 7r ) + 
VKM^M^ + Oia 



3 dM n + 

2\ 



1 



--M-n 



(M w+ - M w o) + 



2 dVl 



3 dM„o 



M t0 ) + O(c 
(4 



i.e. to order 0(a) there is no isospin breaking in the isoscalar TPE. For the isovector TPE, we have 
the general structure 



V, 



2- 



rf rf V 2 \ (M w+ , M n+ ) + (n • r 2 - rf rf) V 2 \ (M w+ , M w o ) 

J V 2 \(M n + , M^+) for pp and nn 

\ 2V 2 \{M^ + ,M^)-V 2 \{M^+,M n+ ) for np,T = 1 



(4.10) 



Clearly, for the (and the nn) system, we have to use the charged pion mass. In case of np, we 
use for the numerical evaluation the expression in eq.( 4,10| ). If one again considers the pion mass 
difference as a small effect and Taylor expands around its average value, one finds that to order a 
one has to use the neutral pion mass in V 2 \ for the T = 1 np system. 



The 7T7 exchange diagrams have already been obtained in ref.jl4| and we use the results obtained in 
that paper omitting all computational details here (although it should be stressed that the treatment 
of the divergent loop integrals has to be considered approximative, however, for the accuracy of our 
calculation it is sufficient). Due to isospin only charged pion exchange can contribute to the n'j 
potential V n ^ and thus it only affects the np system. The potential has the form 



9\ 



(1-/3 2 ) 2 
2/34(1 +/F 



Ty ■ t 2 - r{ t 2 ) ffi ■ qa 2 ■ q V nl ((3) 



ln(l + /3 2 ) + -^ 



1 

2/3 2 



27 



l + P 2 



(4.11) 



Here, (3 = \q\/M 7r + and 7 is a regularization scheme dependent constant. Consistent with the renor- 
malization procedure adopted in ref.Q (no finite renormalization) , we set 7 = and the expression 
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for the potential simplifies accordingly. The analytical form of V n ^ is similar to the one of the OPEP 
but it differs in strength by the factor a/n ~ 1/400. Due to this inherent smallness, one only expects 
it to have some influence on the S-wave scattering length. 

Finally, we consider the contact interactions. In the isospin limit and at leading order (i.e. with 
no derivatives), we have two independent contact interactions and seven at NLO (two derivatives). 
Particular linear combinations of these terms can be projected onto the appropriate partial waves, 
specifically the two S-waves (So, 3 Si), the four P-waves ( 1 -Po ) 3 Pj) and the mixing parameter e\ of 
the coupled 3 Si — 3 D\ system (for T = 0). For the two nucleon systems with T = 1, one has of 
course only the singlet S-wave and the three triplet P-waves. In case of CIB and CSB, each of the 
next-to-leading order S-wave contact terms consists actually of three different contributions, namely 
the isospin symmetric part C sym (with two derivatives) as defined in ref.Q, the CIB part ~ r 3 r| and, 
the CSB part ~ (r 3 + r| ) . Four the P-waves, we only have the two-derivative isospin symmetric 
contact terms since the equivalent CIB and CSB contact interactions only appear at two orders higher. 
Therefore, we concentrate on the S-waves. The evaluation of these interactions in the different isospin 
channels leads to the following contribution to the various NN potentials 

C pp = C sym + Ccib + Ccsb , (4-12) 

Cnn = Csym + CdB ~ CcSB , (4-13) 

C np = C sym — Ccib • (4-14) 

As long as only the pp and the np phase shift are used to pin down the parameters accompanying 
these contact interactions, one can extract from the contact interactions 

Cpp — C np = 2Ccib + Ccsb , 

Cpp + C np = 2C sym + Ccsb • (4-15) 

Only if one has access to a third observable, a separation of these three different contributions is 
possible. In the *So channel, the nn scattering length and the effective range come into play. These 
can be measured by various methods (as discussed below) and for that particular partial wave, we 
can then separate the isospin symmetric, CIB and CSB contributions to the coupling constant C via 

1 

Ccsb = -^{C vv - C nn ) , (4.16) 

Ccib = -^{C PP — C np — Ccsb) • (4-17) 
This finishes our formal discussion and we now turn to the results. 



5 Results and discussion 
5.1 Fitting procedure 

First, we must fix the parameters. Throughout, we work with g& = 1.26, m p = 938.27 MeV, m n = 
939.57 MeV, M n o = 134.98 MeV and M n + = 139.57 MeV. Next, we describe the fitting procedure 
to pin down the LECs. We proceed along the same lines as detailed in Q. We perform two types 
of fits. The global ones are obtained by fitting to the phase shifts of the Nijmegen partial wave 
analysis for the np and the pp case [p6f| . More precisely, the isospin-symmetric LO and NLO LECs 
are determined from a best fit to the S-wave 1 So and the three triplet P-waves 3 -Po,i,2 at laboratory 
energies of 1, 5, 10, 25 and 50 MeV and the NLO LECs entirely from the singlet S-wave (for np 



14 



and pp) .p 5 ! Alternatively, for the S-wave we also use the scattering lengths to determine the isospin 
breaking LECs. This fit is marked by a We work with a sharp momentum-space regulator A 



in the LS equation. The errors for the partial waves are taken from ref. [26|. For a certain range of 
A, the x 2 distribution is shallow as expected in any cut-off EFT. Here, we are only interested in the 
cut-off range which allows to simultaneously fit the np and the pp systems. This is the case for the 
range 300 MeV < A < 500 MeV. The upper limit in this range is given by the fit to the pp partial 
waves. We remark that as a check we have reproduced the results of ref.Q for the np case (for the 
parameters and cut-offs used there). The corresponding values for the LECs are collected in table [l] 
for A = 500 MeV and the cut-off dependence of the LO and NLO S-wave LECs is given in table 0. 





So 




3 Po 


3 Pi 


3 P 2 


C£° [10 4 GeV' 2 ] 
C* LO [10 4 GeV" 4 ] 


-0.133 
1.822 


-0.093 
2.125 


1.335 


-0.394 


0.191 


C p L ° [10 4 GeV' 2 ] 
C^ p LO [10 4 GeV" 4 ] 


-0.129 
1.822 


-0.086 
2.125 


1.335 


-0.394 


0.191 



Table 1: Values of the LECs for in the various partial waves for a sharp cut-off A = 500 MeV. 



A [MeV] 


300 


400 


500 


C n L ° [10 4 GeV' 2 ] 
[10 4 GeV" 4 ] 


-0.1550 
1.613 


-0.1545 
1.621 


-0.1331 
1.822 


C£° [10 4 GeV- 2 ] 
CNLO [i 4 GeV" 4 ] 


-0.1537 
1.613 


-0.1532 
1.621 


-0.1290 
1.822 



Table 2: Values of the LECs for in the So partial waves as a function of the sharp cut-off. 



5.2 Phase shifts 

Having determined the LECs, we can now predict the S- and P-waves for energies larger than 
50 MeV and all other partial waves are parameter free predictions. To be specific, we employ the 
cut-off A = 500 MeV. In Fig||, we show the 1 S partial wave for energies up to 200 MeV. As noted in 
ref. [||], for an improved description of this partial wave up to E\ a b = 300 MeV, one has to go to NNLO. 
The triplet P-waves are collected in Fig.||. As already pointed out in 0], the pion mass difference 
plays an important role in 3 Po for which we get a good description up to E^b = 200 MeV. The good 
description of the 3 Pi phase up to the inelastic threshold can be traced back to the dominance of 
the OPE. The insufficient description of the 3 P2 wave was already be observed in [Q] and can be 
cured to some extend at N 2 LO. Some higher partial waves ( 1 D2, 3 E2, 3 -F3, 1 G4), which are free of 
tunable parameters at NLO, are shown in Fig.[|. Depending on the partial wave, we find satisfactory 
description of these waves up to 150 MeV. This holds in particular for the mixing parameter €2 not 
shown in the figure. These results can be improved by going to NNLO, but for our purpose the 
accuracy achieved for lab energies up to 150. . .200 MeV is sufficient. In particular, it allows to study 
CIB and CSB which is most pronounced in the threshold region. 

#9 We are aware that in the absence of magnetic interactions and vacuum polarisation the first energy in these fits is 
somewhat problematic. 
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5.3 Charge independence breaking 



In the absence of nn data, we define 



Xcm = X pp — X, 



X = {5,a,r,v 1 ,v 2 ,v 3 } 



(5.1) 



for the pertinent phase shift, scattering length and range parameters (in a given partial wave). We 
show first the energy dependence of 5cm in the various partial waves, in comparison to the Nijmegen 
PSA and the most recent results obtained from the CD-Bonn potential. The S- and P-waves are 
collected in Figs.^ and [?], respectively. In the S-wave, our prediction for energies above 50 MeV 
shows an interesting difference to the result of the CD-Bonn potential (filled squares in fig||. In that 
approach a phenomenological term had to be added to describe correctly the CIB in the scattering 
lengths, this affects in particular the CIB difference above I?i a b = 50 MeV as depicted by the open 
squares. The prediction by ignoring this additional term, whose explicit structure is not spelled out 
in ref.0, agrees with our full NL0 result. For 3 -Pi,2) our predictions based on the Nijmegen PSA are 
distinctively different from the CD-Bonn results, however, these phase differences are in both cases of 
very small magnitude. The CIB in the higher partial waves is shown in Fig|| only in 3 F% one finds 
some difference between our prediction and the one obtained using the CD-Bonn potential. 

We now turn to the scattering lengths and range parameters. To be specific, we consider the 1 Sq 
partial wave. The np and pp scattering lengths as a function of the cut-off A are given in table § 



While the result for a 



np is, of course, independent of the cut-off (within an accuracy of about 0.3%), 

It is, 



the pp scattering length varies due to the subtraction ambiguity discussed in section |3.4 
however, gratifying to see that within the range of allowed cut-offs, this variation is modest 
\Sa P p/a P p\ ~ 0.6 fm/17 fm ~ 3%. Consequently, we get 



i.e. 



«CIB = 6.14 . . . 6.75 fm 



(5.2) 



for A = 300 . . . 500 MeV. 



A [MeV] 


300 


400 


500 


a np [fm] 
a pp [fm] 


-23.65 
-17.51 


-23.67 
-17.53 


-23.71 
-16.96 



Table 3: np and pp scattering lengths in 1 So partial wave for the range of cut-off parameters discussed 
in the text. 



To further dissect the physics underlying CIB, we list in table || the various contributions to the CIB 
in the scattering length and effective ranges due to the neutron-proton mass difference (Amjy), the 
pion-mass difference (AM,,-) in OPE and TPE, from nj graphs and due to the contact interactions. 
These contributions have been obtained in the following manner: The total potential is a sum of 
different terms and to assess the importance of one particular mechanism, we simply switch that 
term off and obtain its contribution as the difference between the full result and the one in the 
absence of that term. Due to the inherent non-linearity of the problem (iteration of the potential, 
pseudo-boundstate close to threshold) the sum of the individual contributions does not give the full 
result. In the literature one can find a set of prescriptions to rescale these contributions based on an 
expansion in the inverse of the scattering length. We do not follow these schemes but stick to the 
prescription just outlined. As expected, the pion mass difference is a very important effect, but is 
alone insufficient to account for the CIB in the scattering lengths. The TPE contribution deserves 
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some discussion. It has already been worked out in ref . 1 15 ] , where a similar result with an opposite 
sign was found. There are several reasons for that disagreement. First of all, the authors of ref. [15] 
performed all calculations in configuration space whereas we are working in momentum space. Several 
zero range terms arising by Fourier transformation of the non-polynomial part of eq.(|4.5|) have not 



been taken into account in ref. [15]. Secondly, we are using a sharp momentum cut-off whereas some 
smooth local coordinate space cut-off has been applied in ref.|l5|. Finally, the charge independent 



part of the potential was chosen in [15] to be the one of the Argonne AV18 potential. In this sense, 
our description is more consistent and is entirely based on the potential derived within chiral EFT. 
Therefore, a direct comparison with the numerical results of |ll| should not be made. We find the 
same ordering in the importance of the three types of diagrams, box larger than triangle larger than 
football. However, there is also a strong cut-off dependence, e.g. the box diagram contribution 
to ocib changes form -1.6 fm to -1.1 fm as A varies from 500 to 300 MeV. Such a strong cut-off 
dependence has, in fact, to be expected, since neither the contribution to CIB of the complete CIP 
TPE nor of its non-polynomial part, shown in table ||], correspond to any observable quantity. Only 
the complete CIB contribution to phase shifts can be measured experimentally. The irj contribution 
agrees with the findings in ref. [|l4[ . Finally, we would like to stress the rather sizeable contribution of 
the contact interactions to aciB) which deserve some study within models of CIB. 





Amjv 


AM W OPE 


AM W TPE 


7T7 


contact 


a [fm] 


0.1284 


2.9173 


-1.1376 


-0.4203 


4.5083 


r [fm] 


-0.0019 


-0.0804 


0.0135 


0.0066 


-0.0739 


V\ [fm 3 ] 


-0.0018 


-0.2684 


0.0083 


0.0143 


-0.0589 


V2 [fm 5 ] 


-0.0120 


-2.2246 


0.0541 


0.0893 


-0.4249 


v 3 [fm 7 ] 


-0.0719 


-17.800 


0.2830 


0.5606 


-2.5355 



Table 4: Various contributions to CIB in the 1 Sq partial wave for a sharp cut-off A 
explained in the text. 



500 MeV as 



5.4 The nn system and charge symmetry breaking 

To discuss charge symmetry breaking, one has to know the nn scattering length. With the "standard 
value" of a nn = — 18.9 ±0.4 fm one obtains the value given in eq.( |1.3| ). However, a recent measurement 
performed at Bonn University using nd breakup combined with exact solutions of the three-body 
problem gives a sizeably lower value, a nn = —16.4 ± 0.4 fm ^ TC . For the central value of the Bonn 
result, we get 

acsB = -0.9 fm , (5.3) 



which is smaller in size than the value of eq.( |1.3|) and of opposite sign. As discussed before, we can 
not predict this scattering length due to the presence of a four-nucleon LEC but rather use the value 
of a nn as input to predict the nn 1 Sq phase shift shown in Fig.^ in comparison to the np and pp 
phases. Note also that CSB dies out much faster than CIB, at = 25 MeV the two nn curves 
are identical. It is also interesting to investigate the origin of CSB in the EFT. At NLO, there are 
two effects, namely the nucleon mass difference and the short-distance physics encoded in the leading 
order CSB four-nucleon interaction. While in the case based on eq.(|1.3|), the nucleon mass difference 

* 10 There is a mild scatter on this result depending on the energy and normalization which is, however, within the error 
bar given. 
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effect (0.33 fm) and the contact term (1.27 fm) go in the same direction, for a scattering length as 
in eq,(|5.3]), the short distance physics encoded in the contact interaction (-1.21 fm) has to overcome 
the positive contribution from the nucleon mass difference (0.29 fm). This becomes more transparent 
if we work out the LECs as defined in eqs. fl4.16 , 4, 17| ) at leading order. Using the standard value for 



a nn , we have (in units of 10 4 GeV 2 ), 

Com = 0.0026 , C CS b = 0.0012 , (5.4) 

whereas the value of the scattering length from the Bonn measurement gives 

C cm = 0.0036 , Ccsb = -0.0008 . (5.5) 

We remark that the importance of the contact term contributions to CSB is in agreement with findings 
in the CD-Bonn potential, where CSB is largely driven by TPE with non-nucleonic intermediate 
states. Such contributions are subsumed in the contact interactions of our approach. It is also 
interesting to compare these NLO results with ones obtained in the KSW scheme [11], cf. Fig. 10 



As it has already been found in the isospin symmetric calculations, both schemes give similar results 
in the low energy region but the KSW scheme can not be applied systematically for momenta larger 
than the pion mass. 



6 Summary 

In this paper, we have calculated charge independence and charge symmetry breaking in the two- 
nucleon system based on a chiral effective field theory. This naturally extends the results for the np 
system presented in refs.[||, |j. The results of this investigation can be summarized as follows: 

1) Based on a modified Weinberg power counting (as explained in ref.||), we have systematically 
included strong and electromagnetic isospin violation in a chiral two-nucleon potential at NLO. 
Strong isospin violation is due to the light quark mass difference and can easily be incorporated 
by means of an external scalar source. The electromagnetic effects due to hard photons are 
represented by local contact interactions. Soft photons appear in loop diagrams and also gener- 
ate the long-range Coulomb potential. The corresponding classification scheme for the various 
contributions to CIB and CSB based on the extended power counting at L0 and NL0 is given 



in section 3.2 



2) The resulting potential consists of two distinct pieces, the strong (nuclear) potential including 
isospin violating effects and the Coulomb potential. The nuclear potential consists of one- 
and two-pion exchange graphs (with different pion and nucleon masses) , iT'j exchange diagrams 
and a set of local four-nucleon operators (some of which are isospin symmetric, some depend 
on the quark charges and some on the quark mass difference). Since the nuclear effective 
potential is naturally formulated in momentum space, we use the matching procedure developed 
in ref.[25| to incorporate the correct asymptotical Coulomb states. The necessary regularization 



of the potential is performed on the level of the Lippmann-Schwinger equation, using a sharp 
momentum space cut-off A. 

3) The low-energy constants accompanying the contact interactions and the cut-off A are deter- 
mined by a simultaneous best fit to the S- and P-waves of the Nijmegen phase shift analysis in 
the np and the pp systems for laboratory energies below 50 MeV. This allows to predict these 
partial waves at higher energies and all higher partial waves. Most physical observables come 
out independent of the cut-off for A between 300 and 500 MeV. The upper limit on this range 
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is determined by the pp interactions. The resulting phase shifts are shown in figs||-[|, and the 
CIB phases 5qib = $pp — $np in figs.||-||, in comparison to the outcome of the Nijmegen PSA and 
recent results based on the CD-Bonn potential Q. 

4) We have studied in detail the range expansion for the np and the pp system. For the range of cut- 
offs, the pp scattering length varies modestly with A due to the scheme-dependent separation 
of the nuclear and the Coulomb part, \5a P p/a pp \ ~ 0.6 fm/17 fm ~ 3%. We have dissected the 
various contributions to CIB in the scattering length, reconfirming the importance of the pion 
mass difference in the OPE. The TPE contribution is smaller in size but of opposite sign, and 
also strongly cut-off dependent. 

5) We have also studied the nn system, performing effective range fits based on the "standard" 
value for a nn = — 18.9fm and using the most recent result from deuteron break-up, a nn = 
— 16.4fm M. This leads to rather different physics underlying charge symmetry breaking, as 



discussed in section 5.4 



We have shown that charge independence and charge symmetry breaking can be systematically studied 
in an effective field theory based on a modified Weinberg power counting. In the future, it would be 
of interest to also study the effects of two-pion exchange diagrams with the inclusion of explicit delta 
degrees of freedom to further dissect the physics underlying CIB and CSB. In addition, one should 
also include vacuum polarisation and magnetic moment corrections to the Coulomb potential if one 
wants to achieve the same precision as in the so-called modern potentials. Work along such lines is 
in progress. 
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A Closer look at the two— pion exchange potential 



In this appendix we discuss in more detail the CIB contributions from the TPEP. In fact, we will 
show that one generates some polynomial pieces at NL0 which lead to finite shifts in some LECs. 
However, the effect of these shifts will turn out to be numerically very small. Nevertheless, when 
improving the accuracy of our potential, such effects have to be taken into account. 

In the isospin symmetric case, we can write the NLO TPEP V^\_ lo in the form, see also eq.( f4.5| ), 



Vr 



(2) 



2ir,l— loop 



V, 



TPEP 

non— polynom 



+ (S 1 +S 2 q 2 )(f 1 -T 2 ) + S 3 



(<7i • q) (<r 2 • q) - (<7i • a 2 )q* 



(A.l) 



with 



-.mh^ {- 18M -( 5 9A - 2<A) ^ ^ - Ml{Qlg\ - Ug\ + 4) 

+l8M*(5g 4 A - 2g 2 A )J - 3(3<?i - 2g A ) J 2 



(A.2) 
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52 = ^^{(- 2 35Ul09i + l)ln^-^(13^ + 2<7i) + (23^-10^-l)J }, (A.3) 

5 3 = - ca ot,A ln ~ + o - J » ( A - 4 ) 



T/TPEP _ t/0 , t/1 -r • -f / A c\ 

^polynom 'polynom ~r 'polynom ' 1 ' ' 2 • v-^- - "/ 



64tt 2 F4 1 e 3 

in terms of the divergent loop functions Jq,2, the IR regulator e (for details, see ref.[Q]) and the non- 
polynomial part is given in eq.( |4.5| ). From here on, we concentrate on the polynomial terms, which 
have the genuine structure 

'polynom '"'polynom ''polynom 

Again, the pion mass difference AM% ~ a will lead to CIB. Following the same arguments as used 
in section ||, up to corrections of order a 2 , one uses the various pion masses in the polynomial pieces 
as explained for the non-polynomial ones. To be specific, consider the T = 1 np case. We denote 
the corresponding LECs by Cs,CT,Ci. Putting in the appropriate pion masses, one finds for the 
isovector piece 

^polynom^p, nn ^polynom, np ^polynom (M n ± ) ^polynom ("^vr ) 

= {A + Bq 2 )^ -f 2 . (A.6) 

The term ~ A leads to a relative shift between Cp®, C„„ , and C^® in the 1 So wave. However, we 
do not need to know the precise value of this shift (and thus of the constant A) since these LECs are 
any way fitted independently. To the contrary, one has to know the value of B, 

B ~ 384^ {-™<A + «*1 + l) = CV . (A.7) 

This leads to finite shifts in the NLO LECs, 

a ^yNLO _ ^NLO ,-yNLO ( a o\ 

LAkja L i — ^pp,nn, 3 Lj ~ u np, s Lj > V A -°J 



in the usual spectroscopic notation. We find 



AC^°=4irB, AC 3 N P LO = -^B, (AS 



so that 



= [C^? So - 0.0063] • 10 4 GeV" 4 , 
CgS. 3P, = fe°P, + 0-0042) • 10 4 GeV" 4 . (A.10) 



' pp,nn, 3 Pj I np, Pj 

The effect of these shifts is numerically very small, as demonstrated for the 3 P and 3 P 1 waves 
in fig. 11. 
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Figures 
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Figure 1: Phase shifts for the pure Coulomb potential at E^ = 1, 10, 100, 300 MeV as a function of 
the matching radius R. 
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Figure 2: 1 S'o phase shift for Reid93 plus Coulomb potential at = 10 MeV (solid line) for varying 
matching radius (dashed line). 
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Figure 3: 1 So phase shifts for the np and pp systems (solid and dashed line, respectively) im compar- 
ison to the Nijmegen PSA, filled circles: np, open circles: pp. 
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Figure 4: P-wave phase shifts for the np and pp systems im comparison to the Nijmegen PSA. For 
notations, see Fig.||. 
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Figure 5: Higher partial waves for the np and pp systems in comparison to the Nijmegen PSA. For 
notations, see Fig.||. 
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Figure 6: CIB 1 Sq phase shift (solid line) im comparison to the Nijmegen PSA (dashed line) and the 
CD-Bonn potential (filled and open squares, as explained in the text). 
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Figure 7: CIB P-wave phase shifts in comparison to the Nijmegen PSA and the CD-Bonn potential. 
For notations, see Fig.B. 
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Figure 8: CIB phase shifts in the higher partial waves in comparison to the Nijmegen PSA and the 
CB-Bonn potential. For notations, see Fig.||. 
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Figure 9: Range fit for the nn 1 Sq phase shift based on different scattering lengths (a nn = — 18.9fm: 
long-dashed line, a nn = — 16.4fm: short-dashed line) in comparison to the np and pp phases. For 
further notations, see Fig.||. 




Figure 10: Range fit for the nn and np So phase shifts. The long-dashed and long-short-dashed lines 
refer to the np and nn results at NLO in the KSW counting [jllj] . The results of the present approach 
are: np solid line, nn with a nn = — 18.9, —16.4 fm: dotted and and short-dashed line, respectively. 
The filled circles are the Nijmegen PSA values. 
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Figure 11: Effect of the polynomial TPEP shifts on the ^o, 3 P , 3 Pi waves as discussed in the 
appendix. The (solid) dashed lines correspond to (not) incorporating the finite shifts as given in 
eq. flA.10D . 
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